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Abstract. We perform a Random Matrix Theory (RMT) analysis of the 
quantum four-state chiral Potts chain for different sizes of the chain up to size L = 
8. Our analysis gives clear evidence of a Gaussian Orthogonal Ensemble statistics, 
suggesting the existence of a generalized time-reversal invariance. Furthermore 
a change from the (generic) GOE distribution to a Poisson distribution occurs 
when the integrability conditions are met. The chiral Potts model is known to 
correspond to a (star-triangle) integrability associated with curves of genus higher 
than zero or one. Therefore, the RMT analysis can also be seen as a detector of 
"higher genus integrability". 
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Introduction 

Initially developed in the framework of nuclear physics or atomic physics , Random 
Matrix Theory (RMT) provides a versatile characterization of chaos Q. Since the 
pioneering work of Wigner H Dyson (§, and Mehta [§ RMT has been applied 
successfully to various domains of physics. As a limiting case RMT signals the 
emergence of integrability, which shows up in the change of the generic Wignerian 
level spacing distribution into Poissonian or Dirac distributions. The first examples of 
this appeared when considering simple harmonic oscillators (totally rigid spectrum) 
or free fermions models || [j], ||. The reduction to Poisson distribution reflects 
nothing but the independence of the eigenvalues. This change in the distribution 
may sometimes come from a dimensional reduction of the model, like in the so-called 
disorder solutions ]|, pT| . It can also be found in genuinely correlated systems, 
the reduction being now associated to Bethe Ansatz integrability [H| |l5|] or 
Yang-Baxter integrability [jl6| with rational or elliptic functions. It is natural to ask 
whether this link between Poisson reduction and Yang-Baxter integrability still holds 
when the solutions of the Yang-Baxter equations are no longer parametrized in terms 
of abelian varieties. The perfect example to address this question is the chiral Potts 
model for which Au-Yang et al. found a higher genus solution Jl^] of the Yang-Baxter 
equations. These solutions appeared in the two-dimensional classical chiral Potts 
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model on an anisotropic square lattice (see for instance JT^]). As a consequence of 
the Yang-Baxter equations, there exists a family of commuting transfer matrices, also 
commuting with some quantum hamiltonian given below for the chiral Potts model 
(see (Q) in the following). The RMT approach can be applied directly to analyzing 
the eigenvalues of the transfer matrices of the two dimensional classical models ||, 
but of course it is much simpler when applied to quantum hamiltonians, since the 
latterfj] does not depend on the spectral parameters. It is also numerically much more 
convenient. For the sake of simplicity we restrict ourself to the RMT analysis of the 
quantum hamiltonian. 

The paper is organized as follows: in section [j] we recall some results about the 
chiral Potts model. In section || we review how to use in practice RMT in the context 
of quantum hamiltonian. In section || we review some symmetries of the quantum 
hamiltonian of the chiral Potts model, and discuss time-reversal invariance. Our 
numerical results are presented in section [|, where we also discuss the unexpected 
occurrence of a GOE statistics. 



1. The quantum chiral Potts chain. 

The Hamiltonian of the quantum chiral Potts chain first introduced by Howes, 
Kadanoff and den Nijs and von Gehlen and Rittenberg |2(J is defined by : 

JV-l 

H = H X +H ZZ = J2 H n+i =~EE ^< X o) n +^n{Z j Z} +i r}{l) 

j j n=l 

where Xj = In ® • • • ® X ® ■ ■ ■ ® In and Zj = Ijv® ,,, ®Z'<8>- ,, (8 In- Here In is a 
N x N unit matrix while X and Z are N x N matrices, in j-th position in the tensor 
product, with entries : 

Zj, m = Sj, m exp[2m{j — l)/N] and X^ m — 5j, m +i (mod N) 

The self-dual model |2l[] corresponds to a n = aTT- Conformal theory analysis of 
the 3-state model can be found in p2| . Some spectral analysis of this model have been 
performed for the quantum self-dual model or the 3-state model [ pp| , p3| , pi] . 

In this paper we restrict ourself to the N = 4 (four-state chiral Potts model) non 
self-dual case. The integrability conditions read (see equations (33a), (33b), (33c) and 
(33d) in 10) : 

^ ^ (2) 



ai 2 +Q3 2 _ aj + al ,^ 

02 Oil 

{al - a 2 3 )(2a 2 2 - ai a 3 ) =0 (4) 

(en 2 - a? ) Cm 2 - mat) = (5) 

There are several simple solutions like (up to a multiplicative common factor) : 

(ai, a 2 , a 3 , al, a^, ai) = (r, 1, r, ±r, 1, ±r) (6) 

or: = (r, 1, r, ±i • r, —1, — ± i ■ r) 

f The fact that the Hamiltonian does not depend on the spectral parameters does not necessarily 
mean that it is blind to the abelian or non abelian character of the integrability varieties. This 
appears when one tries to build the eigenvectors of the quantum Hamiltonian of the quantum chiral 
Potts chain, via Bethe Ansatz, since the method fails for higher genus spectral curves. 
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One of these solutions is a self-dual solution and the other ones are also quite trivial. 
In the last two equations (||) and (||) we choose the second factor of the left-hand side 
to be zero. 

In order to have a real spectrum we choose a,\ = ct2 = a\, 02 = S2* 
and ol = 03* (where the star denotes the complex conjugate) yielding a hermitian 
Hamiltonian. A possible parametrization is then : 

ai = a\ = Vl + r + - r , a% = 1 (7) 

Oi = S3* = \/ n 2 + rn + i y n 2 — rn , 012 = n 

where r and n are real and such that \r\ < \n\ and \r\ < 1. The value n — 1 yields 
the self dual situation. Note that we can scale ct\, ct2, Q.3, a!, 02 and 03 by the same 
common factor which enables to normalize a.2 = 1 in (m). 



2. The RMT machinery. 

Performing a RMT analysis means that one considers the spectrum of the quantum 
Hamiltonian, or of the transfer matrix, as a collection of numbers, and looks for some 
possibly universal statistical properties of this collection of numbers. Indeed, neither 
the raw spectrum, nor the raw level spacing distribution, have any universal property. 
In order to uncover them, one has to perform some normalization of the spectrum: the 
so-called unfolding operation. This amounts to making the local density of eigenvalues 
of the spectrum equal to unity everywhere |25, |2^, E8|. In other words, one 



subtracts the regular part from the integrated density of states and considers only the 
fluctuations. It is believed that the unfolded spectra of many quantum systems are 
very close to one of four archetypal situations described by four statistical ensembles 
emerging from the analysis of the (real) spectrum of random|| matrices ||. For 
integrable models this is the statistical ensemble of diagonal random matrices, while 
for non-integrable systems this can be the Gaussian Orthogonal Ensemble (GOE), 
the Gaussian Unitary Ensemble (GUE), or the Gaussian Symplectic Ensemble (GSE), 
depending on the symmetries of the model under consideration. One-dimensional 
quantum systems for which the Bethe ansatz works have a level spacing distribution 
close to a Poissonian (exponential) distribution [p0| , P{ s ) — cx p(~ s ), whereas if the 
Bethe ansatz does not work, the level spacing distribution can be approximated, if the 
hamiltonian is time- reversal invariant, either by the Wigner surmise for the Gaussian 
Orthogonal Ensemble (GOE): 

P GOE( s )- fsexp(-™ 2 /4) (8) 
or by the Gaussian Symplectic Ensemble (GSE): 

P GSE ( S )^i?Vexp(-i? S 2 ) (9) 

where B = (§) \ -2.263. Note that GOE can also occur in a slightly more general 
framework ("false" time-reversal violation, yl-adapted basis J3l)). When one does not 
have any time-reversal symmetry, the Gaussian Unitary Ensemble distribution should 
appear : 

32 

P GUE( S ) - ^2 s2 ex P (-4 S >) (10) 

§ By random matrices one means that the entries of the matrices are independent Gaussian random 
variables. This is a crucial assumption. Of course, if the entries are not independent Gaussian random 
variables, one can get all kinds of cross-overs between these four statistical ensembles psF. 
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The three above expressions are good approximations of the exact P(s), the latter 
being solutions of particular Painleve equations |52], |33|, |J, [35|, |36| . 

Two-dimensional quantum spin systems were numerically shown to yield GOE 
distribution |2j], |37], [38). Other statistical properties may also be studied, like 
correlations between eigenvalues (see section (^J)), but the core of the analysis will be 
to compare the level spacing distribution of the unfolded spectrum with the Poisson 
and the three Gaussian distributions. 



2.1. The unfolding procedure. 

The unfolding can be achieved by different means ^8|. There is however no rigorous 
prescription and the best criterion is the insensitivity of the final result to the method 
employed and/or to the parameters which any unfolding method introduces, for 
reasonable variation. We denote Ei the raw eigenvalues and q the corresponding 
"unfolded" eigenvalues. The unfolding requirement is that the local density of the q's 
is equal to one. One needs to compute an averaged integrated density of states p(E) 
from the actual integrated density of states: 

J-oo ■ 

and then we take t.- L — N~p(Ei). In order to compute p(E) from p(E), several methods 
are possible. One can choose a suitable odd integer 2r + 1 of the order of 9-25 and 
then replace each eigenvalue Ei by a local average: 

-. i+r 

1 2r + 1 ^ 3 y ' 

Then p(E) is approximated by the linear interpolation between the points of 
coordinates (E' i: i). Another method consists in replacing each delta peak in p(E) 
by a Gaussian distribution centered at the location of the peak and with a properly 
chosen mean square deviation. There are two ways to choose this variance: one can set 
the same mean square deviation for every peak, or even better, one chooses a different 
mean square deviation for each peak, so that the number of neighboring peaks inside 
half- width of the Gaussian distribution is kept constant. Another method is to discard 
the low frequency components in a Fourier transform of p{E). A detailed explanation 
and tests of these methods of unfolding are given in |39| . Note that all these methods 
require an adjustment parameter (the number r defining the running average, the 
mean square deviation itself or the number of neighboring peaks inside half width 
for Gaussian unfolding, the cut-off for Fourier unfolding). When this adjustment 
parameter is large, the smoothing becomes too efficient, and the fluctuations are 
washed out. By contrast too small an adjustment parameter gives a totally rigid level 
spacing: the unfolded integrated density of states coincides with the raw integrated 
density of states. Out of the three methods, the moving average unfolding is the fastest 
one, but the Gaussian with adapted mean square deviation gives the best results. 
Notice that extremal eigenvalues are discarded since they are affected by finite size 
effects and this introduces another, slightly less pertinent, adjustment parameter. 



2.2. Quantities characterizing the spectrum 

Once the spectrum has been computed, sorted and unfolded, various statistical 
properties of the spectrum can be investigated. The simplest one, which is also the 
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most significant, and the most universal, is the distribution P(s) of level spacings 
s = ej_|_i — 6i between two consecutive unfolded eigenvalues ej and €j+i. The 
distribution P(s) will be compared to an exponential distribution and to the GOE 
Wigner law (g). Usually, a simple visual inspection is sufficient to recognize the 
presence of level repulsion |l||, the main property for non-integrable models. In order 
to quantify the degree of level repulsion, it is convenient to use a parametrized 
distribution which interpolates between the Poisson law and the GOE Wigner law. 
>From the many possible distributions, we have chosen the Brody distribution : 

P p (s) = c lS P expi-cis^ 1 ) (12) 

with 

r fB + 2\~\ 1+f} 

C2= T \JTl) and c i = ( 1 + / 3 ) c 2 (13) 

This distribution turns out to be convenient since its indefinite integral can be 
expressed with elementary functions. It has been widely used in the literature. For 
(3 = 0, this is a simple exponential for the Poisson ensemble, and for = 1, one 
recovers the Wigner distribution for the GOE. Minimizing the quantity: 

/>oo 

m= / (Pp( S ) - p( S )) 2 d S (i4) 

Jo 

yields a value of (3 which characterize the magnitude of level repulsion of the 
distribution P{s). We have always found 4>((3) small. When —0.1 < (3 < 0.2, the 
distribution is close to a Poisson law, while for 0.5 < (3 < 1.2 the distribution is close 
to the Wigner distribution. 

If a distribution is close to the Wigner distribution (resp. the Poisson law), this 
means that the GOE (resp. the Diagonal Matrices Ensemble) correctly describes the 
unfolded spectrum, but only at the level of neighboring eigenvalues. If one wants to 
go a step further in the description of the spectrum (at a less universal level) , it is of 
interest to compute functions involving higher order correlations as for example the 
spectral rigidity H: 

/ -, r a+E/2 \ 

A 3 (E) = (-mm (N(e) - ae - bf de) , (15) 

\*> a ' b J a -E/2 I a 

where (. . ) a denotes averaging over the whole spectrum. This quantity measures the 
deviation from equal spacing. For a totally rigid spectrum, as that of the harmonic 
oscillator, one has A3 SC (_E) = 1/12, for an integrable (Poissonian) system one has 
A^°\E) — E/15, while for the Gaussian Orthogonal Ensemble one has AQ®^(E) = 
^z(log(E) — 0.0687) + 0(£ , ~ 1 ). It has been found that the spectral rigidity of quantum 

spin systems follows A^°\E) in the integrable case and A^^^(E) in the non- 
integrable case. However, in both cases, even though P(s) is in good agreement with 
RMT, deviations from RMT occur for A^iE) at some system dependent point E* . 
This stems from the fact that the rigidity A^(E) probes correlations beyond nearest 
neighbors in contrast to P(s). 



3. Symmetry analysis. 

Some symmetry properties of the chiral Potts model can be found in the literature Jjofl . 
We briefly sketch and discuss here the symmetries of the chiral Potts Hamiltonian 
which we use in our analysis. 
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3.1. First properties of the Hamiltonian. 



S3*, 0.1 and ct2 



The hermiticity conditions of the Hamiltonian (|1|) are ct\ — a^, oFf 
real. They are compatible with the parametrization ((?]). 

In this work we concentrate on the four-state case, for which the operators X and 
Z read : 

°\ 




Note that 



X 



xz 



( 




V 1 



1 

/ 



and 



/I 





V 





-1 










zx. 



XZ A 



Z A X, 



(16) 



(17) 



Let p be the 4x4 (symmetric) matrix : 
/ll 1 1 \ 



1 



P= 2 



Let R be the spin reversal a 
fl 

R = 



(18) 






V 



1 i -1 -i 
1-11-1 
\ 1 -i -1 i 

related to the Z 4 discrete Fourier transform. Note that p is symmetric and unitary. 

— (7 (mod 4) : 

°\ 
1 



/ 

One verifies immediately that R is an involution (R = R^ 1 ), that R = p 2 , and 
that the conjugation by R permutes the 4x4 matrices X and X 3 on one side, and 
matrix Z and Z 3 on the other side, i.e. : 

R-XR- 1 =X 3 , RZ-R- X =Z Z (19) 

Matrices X 2 and Z 2 are invariant by the conjugation by R. We introduce the 4 L x 4 L 
matrix Kr which is the tensor product of matrix R, L times : 

K B , = R®R®R---R®R (20) 

This matrix is symmetric, real and involutive, and therefore also unitary. One easily 
verifies, from jl9|), that Kr commutes with H. Let U p be the unitary matrix : 

U p = p®p®p---p®p, K R = U p -U t p = U 2 (21) 

One may perform the change of basis associated with this unitary matrix U p and the 
hamiltonian (m) becomes : 



Hzxx — Hz 



N-l 

Hxx = - [Wi 



{Z Y+a n -{x}x J+1 ) n ] (22) 



j 



n=l 



X^ = X 3 . Obviously both hamiltonians (H) and 



since pXp^ 1 — Z and pZp~ x 
( p2| ) have the same spectrum. 

Along these lines one should recall the existence of a duality symmetry (see |po| , El| 
for duality for the classical models) exchanging the operators Xj and ZjZj +1 in (|l|). 
The dual hamiltonian is : 



N-l 

j n=l 



(Xjr + a^-iZ.Z^X] 



(23) 
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and the duality amounts to permuting the a„'s and cE^'s in ([[]). It is also hermitian 
for a.\ — «3, ~a\ — 03*, a2 real and 02 real. If one compares the real spectrum of (|l|) 
and ( p^ ) one finds (this has been checked for L = 3, 4, 5, 6, 7, 8) that they have the 
same real spectrum only on the representations which are the most "symmetric" with 
respect to the color ( (c, e) = (0, e) see below). This is reminiscent of the situation 
encountered in [Q. 

3.2. Representation theory. 

Eigenstates with different quantum numbers are uncorrelated. It is necessary to 
compare only eigenvalues of states having the same quantum numbers. This amounts 
to block-diagonalizing the hamiltonian (see for instance page 1710 in Q), and this 
is an essential requirement of the method. Due to lattice symmetries, as well as 
permutation of colors for the chiral four-state Potts model ([!]) , there exist a collection 
of operators S, acting on the same space as the hamiltonian H, which are independent 
of the parameters cti and al, and commute with H : [H(ai, aj), S] — 0. The 
block-diagonalization is done with the help of the character table of irreducible 
representations of the symmetry group. Details can be found i n p8| , |39j. 

In our case, that is for hamiltonian (|l|) or ( |22| ) or even (|23D, the analysis goes 
as follows. Matrix X is nothing but the shift operator for the color. Introduce, for 
a chain (§]) of L sites, the 4 L x 4 L matrix: S x = X®X®X---X®X, which 
shifts simultaneously all the spins by one. Using (|l^ ) one finds that Sx and the 
hamiltonian (|l]) commute. This operator Sx generates the abelian group Z±. As far 
as lattice symmetries are concerned, we assume periodic boundary conditions. We also 
introduce the lattice shift operator of one lattice spacing Si a tt- Because of the periodic 
boundary conditions Si a tt commutes with hamiltonian ([[]). Similarly Si a tt generates 
the abelian group Zl- Obviously Sx and Si a tt commute and the total symmetry group 
is generically the abelian Zl x Z4 group. Note that, because of their chirality, these 
hamiltonians do not commute with the mirror symmetry which exchanges site n with 
site L + 1 — n. Therefore the space symmetry group is not the dihedral symmetry 
group Dl. However, if some additional conditions on the parameters 
verified, the symmetry group Dl can reappear. 

• The hamiltonian Hx is hermitian iff cFf = 03* and 02 = 02*. The lattice space 
symmetry group of H x is always the dihedral group Dl, and its spin symmetry 
group is generically (i.e. for aT ^ 03) the group Z4, and becomes the dihedral 
group D4 when cFf = S3. 

• The hamiltonian Hzz is hermitian iff a\ = a\ and Ui is real. The lattice space 
symmetry group of Hzz is generically the group Zl, and is the dihedral group 
Dl when ot\ — a 3 . The spin symmetry group of Hzz is always the dihedral 
group D4. 

For generic r and n in equations (0), the total symmetry group is Zl x Z±. We 
always restrict ourselves to hermitian hamiltonians. Consequently the 4 L blocks are 
also hermitian and they have only real eigenvalues. The diagonalization is performed 
using standard methods of linear algebra (contained in the LAPACK library ||43|l ). 
The projector used to block diagonalize the hamiltonian are : 

Pe,c=[j:^ u )^(j2i Ck Sx) (24) 
\n=0 / \k=0 / 
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with cu — exp(27ri/L). This formula specifies the notations used in the rest of the 
paper, the representations being indexed by (e, c) with < e < L and < c < 4. 
Keep in mind that this block diagonalization ( p4| ) is valid for (Q) or (p3|). For (p2|), 
the generator Sx of the Z4 group is replaced by another matrix S, similar to Sx, and 
( p4| ) is modified accordingly. We will denote Pzxx this unitary transformation. 

3.3. Time-reversal invariance and beyond : the origin of GOE statistics. 

Th existence of a time reversal invariance of the Hamiltonian changes the generic GUE 
distribution into another distribution H ||, |^, 

The anti-unitary time-reversal operator can be written as the composition of a 
unitary operator K with the complex conjugation C : 

T = KC (25) 

In the standard case if is a tensor product over all sites of the chain of some spin 
operator (see for instance equations (26.13b) in |Q|). In the following we will have 
to use a more general notion of time-reversal invariance : K will not be necessarily a 
tensor product. We will only impose that the unitary operator if is a constant matrix 
which should not depend on the parameters of H . For instance, for hamiltonian (|l[), 
K must be indepe ndent of the q ra 's and a„'s. 



In appendix ( Appendix A.l ) it is shown that the time- reversal invariance of the 
hamiltonian H implies that K must be either a symmetric or an antisymmetric unitary 
matrix, together with the following relation between the unitary operator K and the 
hermitian hamiltonian H: 

H = K H* ■ R- 1 = K ■ H l ■ R- 1 (26) 

or equivalently : 

K ■ H l = H ■ K (27) 

where H* and H l are, respectively, the complex conjugate, and the transpose, of the 
hermitian hamiltonian H. 

• Consider first the case where if is a symmetric and unitary matrix. Any 
symmetric and unitary K can be written (see for instance page 224 of ref jH|) as 
the product of a unitary operator U and its transpose, namely K = U • £7*, and 
thus the time reversal symmetry equation ( p6|) becomes U ■ U l ■ H l — H ■ U ■ XJ l or 
equivalently: 

{U^H -Uf = U- l H -U (28) 

In other words, U defines a change of basis bringing H into a symmetric form 
H/ s ^ = U^ 1 ■ H ■ U. The Hamiltonian H being hermitian, Hr s -\ is also hermitian 
and is therefore real symmetric. Its level spacing distribution should have a Gaussian 
orthogonal ensemble statistics if this real symmetric matrix is generic enough. 

• Consider now the case where relation (|2?|), is verified with an antisymmetric 
unitary matrix K. The order of the matrix is necessarily even, namely 27V, 
otherwise the matrix is singular. One can then perform a unitary change of basis 
H — > U~ x Ji ■ U where U not only belongs to the TV-dimensional symplectic |fl5| 
group Sp(N), but is quaternion real ^7j. In that case the level spacing distribution 
will have a Gaussian Symplectic Ensemble statistics [|[ ||. 
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The time-reversal symmetry is a particular case of invariance of the hamiltonian 
under the action of an anti-unitary^ operator A. It can be shown that, provided 
the hamiltonian has a so-called "A-adapted basis" (which is the case if A is an 
involution), the system may display a Gaussian orthogonal ensemble rather than 
the GUE, even if the hamiltonian has neither time-reversal invariance nor geometric 
symmetry. 

The form of condition (|27j) does not depend on the representation of the 
hamiltonian. Performing a unitary change of basis : H — > H' = U ■ H ■ U^ 1 , 
one gets from (|27| ) that : 

H' ■ K' = K' ■ (H'f with : K' = U-K-U* (29) 

where K' is still a symmetric unitary matrix. Notice that K does not transform by 
conjugation. 

Let P denote the change of basis which block-diagonalizes the Hamiltonians, and 
a, (3 denotes the indices of the blocks. H may be represented by the i? Q 's and 
Kbiock = P ■ K ■ P l given by its blocks K a ^'s. Condition ( p7| ) becomes : 

H a ■ K a . a = K a:a ■ Hi and H a ■ K a3 = K a</3 ■ (30) 

Remark : It is crucial that the unitary operator K is a constant matrix. 
Introducing the unitary matrix V diagonalizing an hermitian operator H, and A 
the diagonal matrix of real eigenvalues of H, one sees that V ■ H ■ V^ 1 = A = A* = 
V* ■ H* ■ (V- 1 )* = V* ■ H* ■ (V*)- 1 or: 

h ■ y- 1 ■ v* = v- 1 ■ v* ■ h* = y- 1 ■ v* ■ H f (31) 

One thus sees that the matrix V^ 1 ■ V* , which is a symmetric unitary matrix, is a 
solution of (p7|). However this matrix strongly depends on the parameters of H. From 
a statistical ensemble point of view this means that the ensemble of hermitian matrices 
cannot be reduced to the ensemble of real hermitian matrices with an independent 
Gaussian distribution for the entries in each case. The symmetric unitary matrix K 
we are looking for, has to be independent of the a„'s and a„'s parameters. 

4. Results. 

One question addressed in this paper is to decide whether or not the RMT analysis can 
detect "higher genus integrability". One should recall that the quantum hamiltonian 
exhibits genus zero integrability for self-dual case (c^ = 57), or free fermions 
integrability for some algebraic conditions. In order to avoid these simple cases of 
integrability and stick to higher genus integrability, we choose to move, in the oti,al 
parameter space, along a trajectory crossing the integrable variety given by (pi) to (||). 
In order to avoid the self-dual case, we choose n ^ 1 and fix r. From these values of r 
and n we deduce the values of a\ = a\ and en = 03* and 02 using the parametrization 
. The trajectory in the parameter space is obtained by varying a%. In the following 
we will always consider the following trajectories : 

ai = a\ = vl + r + i\l — r , 0,2 = t 

57=03*= V ?J 2 + rn + iy/n 2 — rn , 5J = n (32) 

if In this respect we can also recall the work of von Gehlen ^| where a ^-symmetry (a Lee- Yang 
symmetry at zero magnetic field) survives for non-zero magnetic field as an anti-unitary symmetry 
on a non-hermitian Hamiltonian. 
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where t, r and n are real parameters. 

Integrability on this trajectory appears at the value a% = 1, We concentrate on 
the value of the best /3b ro dy deduced from ( |T^ ) as a function of the parameter t = a 2 . 

We have constructed the quantum Hamiltonian of the four state Potts model 
for various chain sizes, up to eight (L = 8), i.e. matrices of size up to 4 s x 4 8 = 
65536 x 65536. Since the size of the Hilbert space grows as 4 L , it is difficult to go much 
further. The results displayed below show that the size L — 8 is sufficient to answer 
the question we addressed. Using the complex characters and projectors associated 
with the group Zl x Z 4 (see (gjj)) we have performed the block diagonalization of the 
hamiltonian. The dimensions of the 8 x 4 = 32 blocks, are labeled by (e, c) which are 
respectively to the space index in ( |24| ) and the color index in ( |24| ) . The dimensions of 
the 8 x 4 = 32 blocks are gathered in the following table : 





e=0 


e=l 


e=2 


e=3 


e=4 


e=5 


e=6 


e=7 


c=0 


2070 


2032 


2060 


2032 


2066 


2032 


2060 


2032 


c=l 


2048 


2048 


2048 


2048 


2048 


2048 


2048 


2048 


c=2 


2064 


2032 


2064 


2032 


2064 


2032 


2064 


2032 


c=3 


2048 


2048 


2048 


2048 


2048 


2048 


2048 


2048 



When L is a prime integer, the dimensions become simpler : all the blocks 
have the same dimensions d a u ~ (4 L — 4)/4/L except the blocks of maximal 
symmetry with respect to the space group Zl : (e, c) = (0, c) which have dimension 

^(0, c) = 1 + d-all- 

We first found for each of the 32 blocks, that the eigenvalues are not degenerate 
in each block, and, furthermore, these blocks are irreducible. We then performed 
the unfolding in each block independently. All these calculations have been checked 
against full diagonalization for small sizes, as well as for special parameter sets yielding 
a real dihedral symmetry group. The behavior in the various blocks is not significantly 
different. This is not totally surprising since the dimensions d a of the various blocks 
are almost equal to the average dimension d a — 4 L / (4xL) = 2048. Nevertheless the 
statistics is better for larger blocks since the influence of the boundary of the spectrum 
and finite size effects are smaller. To get better statistics we have also averaged the 
results of several blocks for the same quantum chain size L. We moreover compared 
the four unfolding procedures, again getting similar results. We display the results on 
the largest size L = 8 for the best unfolding procedure namely the Gaussian unfolding. 
Figure 1 shows two level spacing distributions P(s), for, respectively, representation 
(0,0) and representation (7,3), for r — 0.78, n — 1.7, and t = 1.5, which corresponds 
to ax = a\ = 1.334 + i 0.469, a 2 = 1.5, 57. = a^* = 2.053 + i 1.250 and Si = 1.7. 
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L=8, R=(0,0) and (7,3), n=1 .7, r=0.78, t=1 .5 




Figure 1: Level spacing distribution vs GOE, GUE, GSE and Poisson. 



This figure clearly shows that the level spacing distribution is close to the GOE 
level spacing distribution. The GUE and GSE level spacing distributions are ruled 
out. Very similar results are obtained for the other blocks and for others values of 
the parameters away from the integrability value = t = 1. We may compare the 
Brody and the GOE distributions at r = 0.5, n = 2.1, and t = 1.5, corresponding to 
ai = a\ = 1.225 + i 0.707, a 2 = t = 1.5, a! = 03* = 2.337 + i 1.833 and a$ = 2.1, for 
the block ( 0,0) . Figure 2 shows the level spacing distribution and the corresponding 
Brody fit (2^2) for the (least square) best value found to be Pbrody — 0.99. On the 
same figure the GOE level spacing distribution is also displayed, both curves are almost 
indistinguishable. 



L=8, R=(0,0) n-2.1 , r-0.5, t=1 .5, p=0.E 




Figure 2: Level spacing distribution vs GOE distribution. 
As described in section (0) we can test how close we are from the GOE statistics, 
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considering quantities like the spectral rigidity A 3 (E). This is depicted on figure 3 
where the spectral rigidity for the same data as in Fig. 2, and compared with the 
spectral rigidity of the GOE together with the rigidity of random diagonal matrices 
(Poisson). 

L-8, R«(0,0) n-2.1, r-0.5, t=1 .5, b=0.99 




Figure 3: Spectral rigidity A3 (£7) for n — 2.1, t = 1 vs spectral rigidity of the GOE. 

The agreement with the GOE rigidity is good up to a value of E ~ 6, which 
means that the correlation involving up to six consecutive eigenvalues are properly 
taken into account by the GOE description. 

Figure 4 and 5 display the level spacing distribution and the spectral rigidity for 
the integrable case r — 0.5, n = 2.1, and t = 1 which corresponds to ct\ = = 
1.225 + i0.707, a2=i=l,al = aj* = 2.337 + il.833 and = 2.1. 
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Figure 4: Level spacing distribution vs Poisson distribution. 



Random Matrix Theory and higher genus integrability: the quantum chiral Potts modell3 

Figure 4 shows the level spacing distribution compared to a Poisson distribution, 
and also compared to the GOE level spacing distribution. The best Brody distribution 
approximation of the data is found to be (i = 0.04 using a least square fit. We have 
obtained very similar results with other values of the parameters n and r, and for all 
the 32 blocks separately, when t is kept equal to the (higher genus) integrability value 
t = 1. This extremely good agreement with a Poisson distribution is confirmed by 
the calculations of the spectral rigidity displayed in figure 5. The RMT analysis can 
therefore be used to detect integrability even when the integrability is not associated 
to abelian curves. In other words the independence of eigenvalues (yielding the Poisson 
distribution) is not a consequence of the abelian character of the algebraic varieties 
occurring in the Yang-Baxter equations. 

L-8, R-(O.O) n=2.1, r=0.5, t=1 




Figure 5: Spectral rigidity A3 (£7) for n = 2.1, r = .5, t = 1 vs Poissonian spectral 

rigidity. 

This extremely good agreement with an independent eigenvalues situation is found 
for t = 1 exactly. When t is slightly different from 1, the distribution is no longer 
Poissonian (as shown by figure 6) in agreement with the fact that the Poissonian 
distribution should appear only at the integrability value t = 1: as soon as t is no 
longer equal to 1, the independence of the eigenvalues is lost, and eigenvalue repulsion 
sets in. This is seen on the behavior of P(s) for small s. However, in the vicinity of 
t = 1, the full distribution is not exactly a Wigner distribution (/?b r0( jy = 1), and 
is an intermediate Brody distribution. We interpret this fact as a finite size effect. 
Figure 6 shows the level spacing distribution for exactly the same parameters as in 
figure 5 (L = 8, r = .5, n = 2.1, for the block R — (0,0)) except parameter t which 
is changed into t — 1.05. 
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Figure 6: Level spacing distribution for n — 2.1, r = .5, t = 1.05 vs GOE distribution 
together with a Brody distribution for f3 = .36. 

The best (least square) fitting parameter Pbrody is Pbrody = 0.36. This 
intermediate value, between and 1, is a consequence of the finite size of the quantum 
chain. One can certainly expect fibrody would tend, in the thermodynamic limit, to 
the GOE value (3brod y = 1- In order to quantify this (finite size) transition from 
integrability to chaos, we calculate the best Brody parameter, as a function of the 
parameter t, keeping r and n constant. Figure 7 displays fibrody, as a function of t, for 
all the representations. 




Figure 7: Best ftbrody parameter as a function of parameter t for all the 32 

representations. 

These results confirm a sharp transition from a GOE distribution to a Poisson 
distribution. In the thermodynamic limit one can expect Pbrody to be equal to the 
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GOE value Pbrody — 1 for every value of the parameter t, except at point t = 1, where 

Pbrody 0- 

To make this change of regime more intuitive, we also show, on figure 8, a window 
on the unfolded spectrum, as a function of parameter t, for L — 7. Only twenty five 
unfolded eigenvalues are represented. One sees clearly, on the unfolded spectrum, the 
level repulsion for t ^ 1 and the level repulsion weakening around t = 1. 




0.5 1 1.5 2 2.5 3 3.5 

Figure 8: Window of the unfolded spectrum as a function of parameter t for seven 
sites. Twenty five unfolded eigenvalues among 4 7 = 16384 are given as function of t. 



4-1. Discussion of the occurrence of GOE. 

The results presented above indicate a clear occurrence of a GOE distribution when 
t 7^ 1 . Numerically the previous results hold for each of the 4 L blocks. This certainly 
requires hamiltonian ([[]) to have additional symmetry properties compared to a generic 
hermitian matrix. In the following we will see that a generalization of the time- 
reversal invariance property for the quantum hamiltonian (|l|), namely condition ( p7| ) 
of subsection (3J3), seems to hold. 

One can look for the matrix K of ( p7| ) in any basis, keeping in mind the particular 
transformation rule (E9h . If we examine form (E2h : 



Hzxx = Hz + Hxx 



H 



M 

zxx 



H 



(as 



X X 



where H 



xx 



is a real symmetric matrix 



H 



xx 



o 3 



(Hz 
with 



H 



XX) 



H 



(as) 
XX 



1, 

2 (Ql 



1 f 



"3) 



and Hj?x ls the antisymmetric matrix: 



H 



(as) 
XX 



Eft* 



t 



■ x)x 



)]-^> 2 -(x^; +1 ) 5 

3 



(35) 
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As a consequence of the hermiticity conditions, in particular ct\ — a^* with ct2 

(s) (s) 

real, Hz is a real diagonal matrix. Thus, since H xx is real symmetric, H zxx is also 
real symmetric. In other words, hamiltonian ( p2| ) is real symmetric when on m = 0, 



in which case K = Kr given by (|20[). It is thus not surprising to see the occurrence 
of a GOE level spacing distribution on hamiltonian (Q) when a\ = 013 is real (namely 
r = 1 with parametrization ([?])). 

When ai m 7^ 0) one looks for a matrix K, independent of the a„'s and a n 's, 
which commutes with H% xx and anticommutes with H xx . 

The existence of K implies that the non zero eigenvalues of H xx appear in 
opposite pairs, which we have checked up to size L = 8, 

4-1.1. Small L cases. • For L = 3, the symmetric unitary matrices K satisfying 
( p7|) are not simple tensor productsQ suggesting that we are not in a strict time- 
reversal invariance framework (see for instance equation (26.15), page 331 in |p^| ). 



For L = 3, the 64 x 64 matrices^ K satisfying ( |27| ) are linear combination of twelve 
quite simple involutive permutation matrices with entries equal to or 1. For periodic 
boundary conditions, none of these linear combinations commute with Si at t, the lattice 
shift operator of one lattice shift spacing. 

This non-trivial form of K is confirmed by its expression in the basis which 
block diagonalizes the Hamiltonian: the off-diagonal blocks K a ^, a ^ /3, of 



Kbiock = P ■ K ■ P* in (|3C|) vanish and one can restrict condition (|27|) to each 
block a — (e, c), namely : 

H a ■ K a , a = K a>a ■ H l a (36) 

The off-diagonal blocks K a ,(3 also vanish for hamiltonian (p2|), the unitary 
transformation ( pi| ) being replaced by Pzxx- 

For L = 3 the symmetry group is Z4 x Z3, and there are 4x3 = 12 blocks 
a = (e, c), labelled in short by an index 0, 1, • • • 11. The blocks Ka^'s can be written 
as : 

where the k a a are simple symmetric unitary matrices with as many entries as possible 
normalized to 1, and where the A Q 's are complex numbers of unit modulus. The block 
matrices k a , a are given in appendix (Appendix A. 2) for hamiltonian (Q) or equivalently 



22]). For instance the block corresponding to the "most symmetric" representation, 

+ Seeking for matrix K = M ® M <g> M satisfying (^?|) when 7^ 0, one finds, from the 

commutation of K with Hz, that M must be a symmetric matrix, and from the anticommutation of 
K w ith H^x t na t the only solution is the null matrix. Of course when o^ m = one gets solution 
(gOJ) taken for L = 3. 

* If K\ and K2 are two unitary solutions of (^), K2 ■ K\ X commutes with the family of H's and 
iCT 1 ■ K\ commutes with the family of H t 's. Thus the set of solutions of ( p7| ) is related to the set 
of matrices commuting with H. 
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namely a = (0, 0) reads for hamiltonian (|l|) as well as for (|2 

1 



ko.o 



(38) 



1 
1 
1 
1 
1 

It is important to note that, up to the multiplic ative (unit modulus) factors A Q 's, the 



blocks k a , a given in appendix ([Appendix A,2j ) are unique. This means that in each 



block a there exists no non trivial symmetry operator commuting with the family (|l|) 
of hamiltonians H. 

[S a , Ha] — 0, Vq„, a^, n = 1,2,3 =>■ S a — Identity (39) 

• For L — 4 we get similar results for Hamiltonian (p2|). The most symmetric 
block (0, 0) yields a 20 x 20 involutive permutation matrix fco,o- It is im portant to 
note that all blocks are quite similar to the ones described in Appendix ( Appendix 
A.2|) and are unique up to multiplicative complex of unit modulus. For L = 4, K is & 
256 x 256 symmetric matrix. If one does not impose the unitarity condition, the set 
of all solutions of ( p7j) reads : 

71=15 

K = ^Pn- A n (40) 

n=0 

where the A n 's are 256 x 256 symmetric matrices whose entries are equal to except 
on at most one entry equal to one for each row or column. In contrast with the L = 3 
case the A„'s are singular matrices (det{A n ) = 0), they are not permutation matrices. 
For certain choice of the parameters p„ one gets, from (|40|), a matrix K which is an 
symmetric real matrix with entries or 1, representing an involutive permutation I\. 

• Similar exact calculations of the blocks of matrix Pzxx ■ Kzxx ■ Pzxx have 
been performed for L = 5 and L — 6. Again one finds that the a ^ (3 off-diagonal 
blocks K a ,j3 vanish and that the 4 L blocks K a a are unique up to multiplicative 
complex of unit modulus. 

As far as the (0, 0) block is concerned one also finds that k Q , are (52 x 52 for 
L = 5 and 178 x 178 for L = 6) simple involutive permutation matrices. 

All these results are detailed on a website |il| where the various blocks k a ^ are 
written for L = 3, L = 4, the blocks such that all the entries are or 1 (but no 
root of unity) are given for L = 5, and L = 6 and furthermore the full 64 x 64 and 
256 x 256 K matrices are written for L = 3 and 4. 

For these values of L (L = 3,4,5,6), the block matrices k a ^ a are remarkable 
matrices with entries 0, or 1, or m-th root of unity (m = 4L). 

4-1-2. Conjecture It is difficult to describe all the 4 L blocks K a ^ a (a ^ (0, 0)). It 
might be easier to describe the 4 L x 4 L matrices K satisfying (27) without imposing 
the unitarity condition in a first step. 

We conjecture that the solutions of (|2?j) are linear combinations of 4 L solutions 
which are involutive permutations /„, in the original basis where X and Z are given 
by ©• 



Random Matrix Theory and higher genus integrability: the quantum chiral Potts modell8 



We conjecture moreover that the block diagonalization of H leads simultaneously 
to a block diagonalization of K into 4 L blocks. The unitarity condition on K 
translates into a simple condition on a multiplicative factor for each block (modulus 
equal to one condition) . The choice of these factors is the only indeterminacy. 



4-1-3. Large on m limit : deformation of a quantized spectrum still yielding GOE. It 
is difficult to find a simple closed expression for matrix K satisfying (pTn , for arbitrary 



L and a, m . We may examine the part H^x °f the Hamiltonian and compare with 



the results of level spacing analysis of Hamiltonian (|l|) or ( |22[ ) for large a 



• Matrix H^x vanishes except on a set of rows and columns where its entries 
are equal to or ±1. Unfortunately the subspace corresponding to these rows and 
columns is not an invariant subspace of Hzxx- ^xx an< ^ ^zxx do not commute. 
Furthermore this subspace becomes "quite large" with increasing chain size L. Let us 
consider the dimension d(L) of this subspace as a function of L. If c?(L)/4 L — > 



when L — > oo one could think that condition ( |27[) "tends to be verified" in the 
thermodynamic limit. In fact this is not the case : the ratio of d(L)/4 L is a monotonic 
increasing function of L. For L running from L = 3 to 12, one gets the following 
values, for successive ratio of d(L)/A L : .375, .406, .647, .663, .677. 

More specifically, the antisymmetric matrix H^x nas the same eigenvalues as 
the diagonal matrix H % ™ z : 

HTz = - ^[ZjZ^-Zjz^] (41) 

j 

that is to say the relative integers 0, ±4, ±8, ±12, • • • ± Am. For instance, for 
L = 12, one gets the eigenvalues ±4, ±8, ±12, ±16, ±20, ±24 respectively 3920928, 
1471932, 268752, 21384, 528, 4 times. 

• We may go back to the level spacing distributions and rigidity calculations 
detailed in section (|J), when ai m is large. This is an interesting situation where the 
spectrum of eigenvalues should be (up to the multiplicative factor oti m ) a deformation 
of a set of integers 0, ±4, ±8, ±12, • • • ±4m. Let us consider again L — 8 but for a 
large enough value of ai m = 200. Figure 9 shows the integrated density of eigenvalues 
for L — 8. It is clear that the eigenvalues are (up to a multiplicative factor a>i m ) 
mainly located around the set of relative integers 0, ±4, ±8. 
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Figure 9: Integrated density of eigenvalues 

Figure 10 shows the corresponding level spacing distribution compared to the 
Poisson distribution, to the GOE distribution, and to the best (least square) fit by a 
Brody distribution. The agreement with a GOE statistics is extremely good since one 
gets flbrody — -93 : 

0=0.93 
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Figure 10: Level spacing distribution vs Poisson, GOE distributions and Brody 

distribution for ftbrody = -93 

Recalling the analysis which yields figures 3 and 5 in section (||) , one also can perform 
rigidity calculations in this strong a.; m limit. Similarly to the results displayed in 
figure 3, the rigidity analysis confirms, for each of the 32 representations, this GOE 
distribution. This is a non trivial limit. This strong ai m limit yields a spectrum which 
is a deformation of a spectrum of relative integers. The unfolding procedure yields a 
level spacing distribution which is still GOE! It does not matter that the eigenvalues 
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are concentrated near integers: what matters is the distribution of eigenvalues around 
these integers which still yields the universal GOE level spacing distribution. 

This very good agreement is a strong indication of the GOE character of the level 
spacing distribution of the hermitian hamiltonian (|l|) in general. 



4-2. Strategy for finding new integrable lattice models. 

One may use RMT analysis to find new integrable lattice models, which is extremely 
difficult analytically especially if they are associated to higher genus solutions of the 
Yang-Baxter equations. 

It has been emphasized |p4| , [H^ , |53| that this type of integrability appears 
when the parameters verify very specific algebraic conditions : these conditions 
express that an infinite set of birational transformations degenerates into a finite 
set ||^, |5^, |||, 57, [58||. We thus have a constructive way to find new possible 
integrability conditions ||59|. However verifying that a particular subvariety of the 
parameter space of the model allows the Yang-Baxter (or the generalized star-triangle) 
relations to be satisfied, remains a very involved analytical task. The RMT analysis 
provides us with a numerically efficient way to verify if these algebraic subvarieties 
yield actual integrability conditions. 

One can show that the general four-state classical two-dimensional chiral Potts 
model has a canonical elliptic parameterization. >From this parameterization one 
may write down explicitly the equations of these algebraic varieties, which are the only 
possible locations for the higher genus integrability conditions |50||. Various analysis, 
similar to the one summarized on Figure 7, of Pbrody as a function of parameter t, 
on various trajectories (0) in the parameter space of the quantum hamiltonian (|lj), 
indicate that the integrable variety @j is the only one with higher genus. Of course 
one cannot exclude the existence of higher codimension integrable varieties avoiding 
the trajectories (Q) we have considered. 



5. Conclusion. 



We have performed a RMT analysis of the quantum four-state Potts chain for different 
sizes of the quantum chain, and for different unfolding methods. Our calculations 
unambiguously exhibit a GOE statistics and excludes GUE (and GSE) statistics. 

Our results indicate that there exists, for arbitrary size L, a symmetric unitary 
matrix K, such that K ■ H* = H ■ K. This can be checked exactly for small lattice 
sizes L — 2,3,4,5,6. We conjecture that such a relation exists for all sizes of the 
chain, and for each of the 4L blocks (jslf). The existence of K would account for the 
statistics we find (GOE rather GUE). 

When the hamiltonian becomes integrable our analysis shows the change from 
the (generic) GOE distribution to a Poisson distribution and this reduction does not 
require the spectral curve to be of genus or 1. 

It is thus interesting to combine this RMT approach with more algebraic methods 
developed in previous publications (2^, . These methods will give the algebraic 

subvarieties where a "higher genus" integrability may appear, if any ... (the infinite set 
of these algebraic subvarieties can be obtained exactly for the four-state chiral Potts 
model @|). As we have shown, the change in the level spacing statistics will signal 
integrability, bypassing the difficulties of the analytical approach. 
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Appendix A. Appendix 

Appendix A.l. Generalized time-reversal invariance. 

The anti-unitary time- reversal operator T can be expressed as the product of a unitary 
operator K and the conjugation operator C, namely T — K-C, where T is projectively 
an involution, namely T 2 = A. Id, and where Id denotes the identity operator. The 
factor A being equal to ±1 as a consequence of the unitarity of K. The time-reversal 
operator T must change the time evolution operator according to: 

Te -iHt T -i = e +iHt (A.l) 

or equivalently 

Ter lHt T = A- e +lHt (A.2) 



Expanding (|A.l| ), or (|A.2j ), in the time variable t, one gets for every order n: 
K (H*) n K* = X-H n Vn 
yielding only two equations: 

KK* = A (A.3) 
KH*K* = A • H (A.4) 



For an hermitian hamiltonian, (A.4) becomes using (A.3) : 

H = K-H l -R- 1 (A. 5) 

where H { is the transpose of H. Since the operator if is a unitary one, (A.3) yields : 



K = A. (iT)" 1 = X.K 1 (A.6) 



where K l denotes the transpose of K. Recalling that A = ±1, we see, from (A.6), 
that K must be either a symmetric or an antisymmetric unitary matrix. 



Appendix A.2. Matrix K as block diagonal matrices for L — 3 for hamiltonian ([£ 

Let lj be the third root of unity : lu = —1/2 — iyo/2, and the consider L — 3 case 
with a = 0, 1, 2, • • • , 11 indexing the twelve blocks defined in (pTj). 

For hamiltonian ( |2^ ) and Kbiock = Pzxx ■ Kzxx ■ Pzxx m ( ftCQ , the off-diagonal 
blocks of Kbiock vanish and one finds the following expressions for the diagonal blocks 
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and fcg g = £>3,3, the block A^o is the same as the one for hamiltonian (|l|). 
Furthermore &2,2 is equal to block k\,i where uj is changed into ui 2 , i.e. £2,2(^0 = 
fci ! i(w 2 ). Similarly one gets k^^(uj) = k^{uj 2 ) as well as ks,s(uj) — kfj(uj 2 ) and 
ku,n(uj) = fci ,io(w 2 ). 
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